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ABSTRACT 

For the past few years, we have observed the central half parsec of our Galaxy in the mid- 
infrared from 2.8 to 5. lyum. Our aim is to improve our understanding of the direct environment 
of SgrA^, the supermassive blackhole at the centre of the Milky Way. This work is described 
in the present paper and by Moultaka et al. 2015 (submitted). 

Here, we focus on the study of the spatial distribution of the ice and gas-phase absorp¬ 
tions. We observed the central half parsec with ISAAC spectrograph located at the UT3/VLT 
ESO telescope in Chile. The slit was placed along 22 positions arranged parallel to each other 
to map the region. We built the first data cube in this wavelength range covering the central 
half parsec. The wavelength interval of the used M-band filter ranges from 4.6 to S.ljim. It 
hosts the P- and R- branches of the ro-vibrational transitions of the gaseous ^^CO and ^^CO, 
as well as the absorption band attributed to the ^^CO ice at 4.675yum. Using two calibrators, 
we could disentangle the local from the line-of-sight absorptions and provide a first-order 
estimate of the foreground extinction. We find residual ices and gase-phase CO that can be 
attributed to local absorptions due to material from the interstellar and/or the circumstellar 
medium of the central parsec. Our finding implies temperatures of the order of 10 to 60K 
which is in agreement with the presence of water ices in the region highlighted by Moultaka 
et al. (2004, 2005). 
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1 INTRODUCTION 

The super-massive black hole (SMBH) Sagittarius A* (Sgr A*) at 
the centre of the Milky Way (MW) has a mass of roughly 4 million 
solar masses and is located at a distance of about 8 kpc (e.g., Ghez 
et al. 2008; Gillessen et al. 2009). It is surrounded by a dense and 
massive star cluster (the NSC, Nuclear Star Cluster) which has an 
estimated mass of 3+1.5x10^ solar masses (Launhardt et al. 2002). 
The mid-infrared images from the Spitzer Space Telescope show 
how the NSC stands out as a separate structure at the centre of 
the MW (Stolovy et al. 2006). In the mid-infrared, the central par¬ 
sec shows a spiral structure, called the minispiral, made of dust 
and ionized gas. It is surrounded by a clumpy, circumnuclear ring 
of dense gas and dust, the CND (Glisten et al. 1987). Due to its 
proximity the Galactic Centre can be studied in detail by the di¬ 
rect observation of individual stars, gas and dust. It represents an 
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ideal and unique case allowing one to analyse thoroughly the direct 
environment of a central SMBH. 

The Galactic Centre is obscured mainly by extinction from the 
diffuse interstellar medium (ISM) present along the line-of-sight. 
About less than one third of this extinction arises from the dense 
foreground ISM (Whittet et al. 1997). It is shown that the mean 
value of the visual extinction towards prominent sources within 
the central stellar cluster is about ~ 27mag (e.g. Schodel et al. 
2010, Scoville et al. 2003, Lebofsky, Rieke, Tokunaga, 1982) with 
a rather smooth distribution across the central 10 to 20 arcsec at an 
angular resolution of 2 arcsec (Schodel et al. 2010, see also Scoville 
et al. 2003). 

Nishiyama & Schodel (2013) find that young, massive star 
candidates can be detected throughout the nuclear star cluster. 
Moreover, Eckart et al. (2013, 2014) showed that there is a sub¬ 
stantial number of faint compact infrared excess sources in the cen¬ 
tral stellar cluster. So far, we can only speculate on their nature 
and origin (see also Muzic et al. 2010, Jalali et al. 2013). Many 
may be associated with young stars or proto-stellar aggregates. In 
particular, there is a small but dense cluster of comoving sources 
(IRS13N) located ~3 arcsec west of SgrA* just 0.5 arcsec north 
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of the bright IRS13E cluster of Wolf-Rayet and 0-type stars. They 
are either dust-embedded stars older than a few Myr, or very young, 
dusty stars with ages younger than 1 Myr (Eckart et al. 2004, 2013, 
2014, Muzic et al. 2008). Eckart et al. (2014) present a first Ri¬ 
band identification and proper motions of the IRS13N members, 
the high-velocity dusty S-cluster object (DSO; Eckart et al. 2014, 
also referred to as G2 by Gillessen et al. 2012, 2013ab), and other 
infrared excess sources in the central field. Eor the IRS13N sources 
the common K- and L-band proper motions indicate that they are 
not only broad-band spectroscopically but also dynamically young 
(Muzic et al. 2008). Eor the DSO the colour information indicates 
that it may be a dust-embedded star rather than a cloud. 

In Moultaka et al. (2004, 2005), we studied the distribution of 
water ices and hydrocarbons towards a number of bright sources 
in the centre of the MW. We estimated the line-of-sight extinction 
in the L-band of the spectrum using a novel approach described 
by Moultaka et al. (2004) and we derived intrinsic spectra of the 
sources. The results show that a substantial amount of water ices 
and hydrocarbons are present in the central parsec around the cen¬ 
tral SMBH. 

In Moultaka et al. (2009) we investigated the circumstellar ma¬ 
terial around the brightest dust-enshrouded sources of the central 
stellar cluster using slit-spectroscopy in the spectral range of the 
M filter from 4.6 to 5.1/im. Eor the 15 bright sources studied in 
that paper, the observations resulted in M-band spectra showing 
the vibration-rotational P and R branch absorptions of the gaseous 
(4.666yum) and {A.lljim). The broad P and R branch 
envelopes are separated by a gap at 4.666yum (the v = 1 - 0 band 
centre; Allamandola 1984). In addition to that, we found a strong 
absorption, centreed at 4.675yum, attributed to a mixture of polar 
and apolar CO ices. Applying a similar method to the one by Moul¬ 
taka et al. (2004) to the data presented by Moultaka et al. (2009) 
we performed a first-order correction of the line-of-sight absorp¬ 
tion due to CO-ice and gas. We found residual absorptions of 
the solid- and gas-phases of the CO that can be attributed to local 
material in the minispiral and the circumstellar medium. In combi¬ 
nation with other data we obtained gas masses of the circum-stellar 
shells of the order of 10 “^ and lO^^M©. 

Previously, Moneti, Cernicharo & Pardo (2001) associated the 
strong gaseous CO absorption to material with a bulk kinetic tem¬ 
perature of ~ lOK towards SgrA*. Eabry-Perot measurements by 
Geballe, Baas & Wade (1989) of the CO R(2) and R(5) lines to¬ 
wards IRS 1, 2, 3, 5, 6 , 7 and 8 revealed that most of the absorption 
takes place in the velocity interval between 0 and 75 km/s. The bulk 
of the foreground absorption is most likely taking place at 20 km/s 
and 50 km/s associated with two giant molecular clouds close to 
the Galactic Centre. The CO gas-phase absorption has also been 
observed in IRS 3, IRS 7 and IRS 12 by Geballe (1986) and Mc- 
Eadzean et al. (1989), but the P and R branches were not resolved 
in these low resolution spectra. While emission from the 4.6/im vi¬ 
brational transitions of CO may originate from the inner few au 
of young stars, absorption is a potentially powerful probe of the 
wind/envelope structure at larger distances from embedded young 
stars. Herczeg et al. (2011) detect emission in the 4.6yum CO funda¬ 
mental line from 14 of the YSOs in a sample of 18 objects originat¬ 
ing from the inner discs for the lower luminosity objects and from 
slow outfiows for objects with high bolometric luminosity. 

In contrast to the gaseous CO absorption, the solid CO has 
a single absorption band at 4.675yum (~0.01yum EWHM) that is 
widened and shifted due to ice impurities (see Hagen, Allaman¬ 
dola & Greenberg 1979, 1980). Initial investigation of the solid- 
phase 4.675yum CO absorption towards IRS 12 by McEadzean et al. 


(1989) suggested that it is dominated by foreground material with 
temperature less than 100 K and possibly as cold as only 20 K. 
Potentially young stars that are still embedded in dense molecular 
envelopes are difficult to identify and study because emission from 
a disc can be confused with possible outflow emission or a larger 
envelope that may e.g. be associated with the minispiral. 

In Moultaka et al. (2009), we also detect a broad absorption 
feature at 4.62yum, discovered by Soifer et al. (1978) and resolved 
by Lacy et al. (1984), that is strongest towards the M2 supergiant 
IRS7. It can be associated with the XCN feature due to cold 
grain mantles that contain molecules with C-N bonds resulting 
from UV-photolysis. A similar feature is produced by laboratory 
UV-photolysis of CO and NH 3 (Lacy et al. 1984, Bernstein et al. 
1995) and is assigned to OCN". Van Broekhuizen et al. (2005) 
report on the detection of this feature in the M-band spectra of 
34 deeply embedded young stellar objects (YSOs), observed with 
high signal-to-noise ratio (S/N) and high spectral resolution. Their 
data enabled us to do the first studies of the solid OCN" abundance 
towards a significant number of low-mass YSOs. The authors show 
that the OCN" abundances are presumably due to photochemical 
and surface chemistry formation mechanism. 

In this paper, we present the analysis of a complete imaging 
spectroscopy data cube obtained via slit spectroscopy and covering 
the 4.6-5. lyum M-band spectrum for the inner 0.5 pc of the central 
clusteiQ. In the following section we describe the observations and 
the data reduction. In section 3, we explain how we derive the line- 
of-sight extinction from gaseous and solid CO. In section 4, we 
analyse the spatial distribution of the CO absorption in its gas- and 
solid-phase across the entire central stellar cluster and show that we 
are able to distinguish between the amount of source intrinsic and 
line-of-sight absorption. Einally, in the last section we discuss our 
results and provide a conclusion. 


2 OBSERVATIONS AND DATA REDUCTION 

We used ISAAC spectrograph located at the ESO UT3-VLT tele¬ 
scope to map the central half parsec of our Galaxy in the M-band 
spectroscopic domain. We needed three periods of observations 
(083.C - June 2009, 085.C - June 2010 and 087.C - August 2011) 
to complete our programme. This ended up with the first data cube 
of the region in the M-band (from 4.4 to 5.1/im). In Eigs. [T]and[ 2 | 
we show the integrated and smoothed integrated map, respectively, 
covering the entire spectral range. 

The optical seeing varied from 0.39 to 1.9 arcsec during 
period 083.C, from 0.4 to 0.8 arcsec during period 085.C and from 
0.7 to 1.5 arcsec during period 087.C. We used a slit width of 0.6 
arcsec providing a spectral resolution of ~800 (i.e. corresponding 
to a velocity resolution of Av = 375km/s). 

To map the central half parsec, we needed 22 slit positions 
that we placed parallel to each others. But because of technical 
problems, we were not able to observe the regions shown in 
black in Eig. [3] Therefore, we interpolated the spectra at these 
slit positions to fill in the whole data cube (see Eigs. [T] and [2]). 
Besides, all slits were not observed with the same integration time. 


^ Resulting from ESO VLT observations of program ID numbers 083.C- 
0675A, 085.C-0469A and 087.C-0214A. 
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therefore, the S/N ratio was different from one slit to another. 

In the following, we explain the data reduction and the build¬ 
ing of the data cube, step by step. There are five different steps. 

• First, all the array images were fiatfielded, corrected for cos¬ 

mic rays and for distortions along the axis of dispersion. The sky 
emission was removed using the chopping technique (offered for 
ISAAC) combined with telescope nodding. The chopper throw dis¬ 
tance was 20 arcsec along the slit. Each chopped frame is com- g 

posed of a positive trace image and a negative one (see an example g 

in Fig. lAlb . Telescope nodding results in two consecutive chopped ^ 
frames A and B shifted 20 arcsec relative to each other. The shift Q 
distance is equal to the chopper throw, such that the positive trace 
image in the first image A is located at the same location on the ar¬ 
ray as the negative one in image B (see example in Fig. lAll) . Thus, 
subtracting one chopped frame from the second (A-B) results in an 
image with two negative traces and one positive trace with twice 
the intensity of the negative ones (see Fig. IaB . The goal of this 
subtraction is to increase the Signal-to-Noise ratio and to remove 
the sky emission lines from the spectra. 

• In the second step, we extracted spectra of the bright sources 
from all the frames and wavelength calibrated them using a Xenon- 
Argon lamp and the third order grating. Then, we made a relative 
flux calibration of the spectra and corrected them for telluric lines 
using AOV-type standard stars and the data reduction commands 
from IRAF software. The procedure to remove telluric lines is to 
divide the extracted science spectrum by a modified telluric stan¬ 
dard spectrum. The modified spectrum is the standard star spectrum 
shifted in wavelength to correct for possible errors in the zeropoint 
dispersion and scaled in intensity (following an exponential law 
that involves a scaling parameter to correct for airmass deviations). 

The optimal ’’shift” and ’’scale” parameters (x/ and S/ in the illus¬ 
tration of Fig. IA21 respectively) obtained during the telluric lines 
correction are stored for each of the extracted spectra of a single 
frame. Telluric standard stars were observed as close to the cor¬ 
responding airmass at which the Galactic Centre was observed as 
possible. Hence we minimized the possibility of calibration uncer¬ 
tainties due to the fact that the sky lines can be narrower than our 
spectral resolution elements. 

• In the third step, for each science frame, the mean value of 
the Xi and Si parameters is calculated on all extracted spectra (v/ 
and Si in Fig. lASI) . This allows us to build a modified telluric star 
spectrum shifted and scaled with the optimal parameters allowing 
us to correct all the spectra of a single science frame at a time. The 
resulting shifted and scaled star spectrum is then stacked along the 
slit axis 1024 times, corresponding to the number of detector pixels 
(see example in Fig. lASI) . In the following, we will call the resulting 
image a ’stacked standard star’ frame. 

• The fourth step is to divide each of the science frames flat- 
fielded, corrected for cosmic rays, distorsions and sky lines, by its 
corresponding ’stacked standard star’ frame. This allows an op¬ 
timal correction for telluric lines of all spectra in the totality of 
the science frames without extracting the spectra (see example in 
Fig.lH. The reduced science frames corresponding to the same 
slit position were then added to each other to improve the S/N. This 
results in one final reduced frame per slit position. 

• In the final step, we constructed the data cube using DPUSER 
softwar^ by positioning properly the resulting 22 slit frames (see 

^ Developped by Thomas Ott http://www.mpe.mpg.de/~ott/dpuser/ see 
also Eckart & Duhoux (1991). 



Right Ascension 


Figure 1. M-band integrated map obtained from the observed data cube. 

Fig.l^. To this end, the different slit frames were shifted ade¬ 
quately along the slit axis (i.e. along the declination axis) to recover 
the observed field with the right relative positions of the bright 
sources. The right ascension axis is obtained by positioning the 22 
slit frames correctly next to each others. This step was done by eye 
but the resulting field of the data cube is recovered very success¬ 
fully as one can see in Fig.|4| This figure shows the smoothed map 
of the integrated intensities along the M-band with contours of an 
M-band image obtained independently with ISAAC imager. This 
figure shows that the positions of the sources are well determined 
within a <0.5 arcsec distance accuracy. 

Most of the sky emission was removed in step one, after the 
A-B subtraction between the two consecutive chopped images. 
But in six slit positions, it was very difficult to remove the sky 
properly. This results in the vertical stripes (seen in Figs. [T] and [2]). 
The slit positions are shown in white in Fig. [3] At these positions, 
all the results described in the paper are obtained by interpolating 
the measurements between the neighbouring slits. 

3 THE LINE-OF-SIGHT EXTINCTIONS 

3.1 Correcting the line-of-sight solid ice ^^CO band and dust 
extinction 

In Moultaka et al. (2009) we formerly corrected the M-band spectra 
of 15 bright sources in the central parsec for the foreground contri¬ 
bution of the solid-phase ^^CO absorption. In that paper, we found 
that this absorption, located at 4.675//m, is very prominent in the 
spectrum of IRS 2L. The red source IRS 2L is located in the cen¬ 
tral parsec but outside the minispiral and the CND molecular disk 
(see its location in Fig|4|). Moreover, it is an early-type star (B-type; 
Clenet et al. 2001) showing no CO bandheads in its K-band spec¬ 
trum. We could, therefore, assume that this star is not affected by 
the local absorption and that the ^^CO absorption observed in its 
M-band spectrum is mainly due to the line-of-sight extinction. We 
smoothed the spectrum of this star to derive a shape of the ^^CO 
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Figure 2. M-band integrated map smoothed with a boxcar of radius 0.3 
arcsec and then using a Gaussian of 0.3 arcsec FWHM. 
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Figure 3. M-band integrated map with the locations of the slit positions 
contaminated by the sky emission lines (shown in white stripes) and the 
locations where no data were obtained (shown in black stripes) because of 
bad weather conditions. 

ice absorption free from all other features like the CO gaseous line 
complex and the hydrogen Pf^ emission line at 4.65/im. Then, we 
normalized to one the resulting spectrum at the wavelengths where 
the continuum is expected not to be absorbed. This spectrum repre¬ 
sents the shape of the foreground ice absorption continuum normal¬ 
ized to unity. Let us call it the ’template extinction spectrum of the 
solid CO band’. It is shown in Fig. 6a of Moultaka et al. (2009) and 
in Fig. IA6I This spectrum shows a broad short wavelength shoul¬ 
der due to a possible XCN absorption present in the spectrum of 
IRS 2L (see Fig. 6(a) in Moultaka et al. 2009). Therefore, we can 
consider that the template spectrum also calibrates the foreground 
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Right Ascension 

Figure 4. M-band integrated map obtained with the data cube corrected for 
the line-of-sight extinction due to dust and to the solid-phase absorp¬ 
tion. Contours: M-band image obtained with ISAAC imager of the UT3 
ESO/VLT telescope. This image shows that our reconstructed field is very 
successful where the positions of the sources are determined within a <1 
arcsec distance accuracy. 

XCN absorption feature. 

Here, we assume that the dominant contribution to the extinction 
by the cold foreground material is in fact due to two components: 
1) dust extinction and 2), the ^^CO ice absorption. The first compo¬ 
nent (dust extinction) can be approximated by a constant continuum 
over the wavelength range (we call it k in equations [T] and [3]). This 
component cannot be estimated from spectroscopy (or from our 
data), but its estimate is not necessary for the results presented in 
this paper. Concerning the second component, the ^^CO ice absorp¬ 
tion, its amount can be variable across the region. This is why we 
diluted the template spectrum by an additive constant continuum d 
to allow for a possible variation in the amount of absorption. To de¬ 
termine this constant d, we adopted the following criterion: when 
we divide the observed spectrum of a Galactic Centre source by this 
diluted template, the resulting spectrum is free from the solid CO 
absorption; this means that the resulting spectrum should show a 
non-absorbed continuum. In this case, the gaseous ^^CO and ^^CO 
line complexes should approximate essential features of the theo¬ 
retical spectrum by Moneti et al. (2001) shown by Moultaka et al. 
(2009) at our spectral resolution (the spectrum is also shown here 
in Fig jATl ). This means that the flux in the wavelength range around 
4.665yL/m, between the P and R branches of the ^^CO line complex, 
is at the same level as that at 4J3fim and 4.8yum (see Fig jAVl ). 

In Moultaka et al. (2009), we applied this procedure and divided 
the observed Galactic Centre spectra by a diluted template spec¬ 
trum with the most appropriate constant d for each case, satisfy¬ 
ing the previous criterion. The corrected spectra are shown in Fig.9 
by Moultaka et al. (20090 We found a mean value for the con- 

^ The fluxes of the corrected spectra for foreground extinction, shown in 
Fig. 9 of Moultaka et al. (2009), are smaller than the non corrected spectra 
of figs 2 to 5 of that paper. This should not be the case since the intrinsic 
corrected spectra are not obscured and therefore are brighter. The reason 
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stant diluting continuum J = 3.7 and a median value of 4. The 
Galactic Centre sources studied in that paper are spread all over 
the 0.5pc central region. This suggests that the template spectrum 
diluted with an additive continuum of 4 is a good approximation 
to describe the foreground overall extinction due to the solid 
absorption. 

This motivated us to use a median value of 4 for the diluting con¬ 
stant d to correct our present data cube for this foreground absorp¬ 
tion. It ends up that, if we call ImtrA^ the intrinsic flux, at wavelength 
T, of a given source in the Galactic Centre (or at a given spatial 
pixel in the field of our data cube), lobs a the observed fiux, at wave¬ 
length T, of the source (or at the given pixel) and Euco band ext a^ the 
template extinction spectrum of the solid ^^CO band, then one can 
write: 

hntrA * [4 -|- Eu^Qf^ond ext a] ^ k — lobs A (1) 

where k accounts for the dust extinction continuum. Let us call 
EucodiiA’ the diluted template extinction spectrum of the ice ^^CO 
band : 


E^^COdil A ^ El2(jQbandext A 

(2) 

Then, equation ([T]) becomes 


j ^obsA 

^intr A , 

Ei^codiiA * k 

(3) 


Here, we derive a corrected data cube for the foreground ex¬ 
tinction by dividing the whole cube with the diluted template ex¬ 
tinction spectrum. The data cube is not corrected for the absolute 
extinction due to obscuration by dust (i.e. k - 1). 

All these approximations can be done since the line-of-sight ex¬ 
tinction across the central 0.5 pc is shown to be constant within 
the uncertainties of <0.3 mag (Schodel et al. 2010; Fig.7) i.e. 
AM<0.015 mag (Viehmann et al. 2005). In addition, with Moultaka 
et al. (2009), we found only a little variation of the additive con¬ 
stant continuum across the region; moreover, the corrected spectra 
of the bright sources obtained in that paper (shown in their Fig. 9) 
are consistently different for different source nature. In particular, 
the corrected spectra of the Helium stars outside the minispiral area 
show no ice absorption (e.g. IRS 16NE) while the dust embedded 
sources in the minispiral have a prominent ice absorption in their 
spectra (e.g. IRS 21). If the residual ice absorption in the corrected 
spectra were due to foreground clouds, one would expect that it is 
also present in the corrected spectra of the non-embedded sources 
like the IRS 16 sources and not only in the dust embedded ones. 

3.2 The foreground R(0) gaseous line correction 

The observed spectra of the bright sources in the central parsec 
show two absorption line complexes (see Moultaka et al. 2009). 
These correspond to the P and R branches of the gaseous 
and rotation-vibration lines. These molecules are probably 
located in molecular clouds along the line-of-sight and in the local 
medium of the Galactic Centre. 

To estimate the amount of these gases in the foreground material 
and in the local medium, we choose to use the optically thin 
isotopic lines for which the optical depth is significantly smaller 

for the small fluxes of fig. 9 is that the absolute value of dust extinction in 
the present wavelength range is not known and cannot be derived from our 
observations. This is a multiplicative factor, called k in the present paper, 
that would rescale the corrected spectra. 


compared to the lines, since the gas is less abundant in 
the ISM (e.g. Geballe et al. 1989, Moneti et al. 2001). In particular, 
we choose the R(0) line located at 4.765 fim where the 

spectra are well corrected for telluric lines and not contaminated 
by sky emission. 

We constructed the optical depth map of the gaseous 
R(0) line using the corrected data cube for the solid 
foreground absorption and dust extinction. The line is not af¬ 
fected by any residual solid CO absorption that is located at bluer 
wavelengths. To derive the optical depths, we approximated the 
continuum by a straight line connecting the spectra from 4.58//m 
to 4.8yum. The optical depth is then derived from the formula: 
^i3cc> = where hontatA and hneatA are the fluxes of the 

continuum and the line at the R(0) absorption wavelength, re¬ 
spectively. 

The optical depths were derived at spatial pixels of the data cube 
where the S/N in the integrated map was higher than 7 on average 
(since each slit position was observed with a different integration 
time). For the remaining pixels, we adopted a zero value. 

The resulting map is shown in Fig. |5] In this map, we interpolated 
the values of the optical depths at the spatial pixels of the slit posi¬ 
tions contaminated by the sky emission lines (i.e. the slit positions 
shown in white in Fig. [3]). 

The resulting map agrees well with our previous work 
described by Moultaka et al. (2009). Indeed, here we find roughly 
the same values of the optical depths observed towards the bright 
sources discussed in that paper. In particular, for IRS 16C, we find 
an optical depth of about 0.10. This value is equal to the one we 
found in our previous paper that was of ~ 0.11 + 0.02. As explained 
in that paper, we use this value to estimate the line-of-sight gaseous 
absorption. Indeed, IRS 16C is an early-type hot Helium star (e.g. 
Krabbe et al. 1995, Najarro et al. 1997) located off the minispiral 
area (see its location in FiglU. Its spectrum shows a pronounced 
R(0) line in a spectral region well corrected for telluric lines. 
Therefore, we could safely assume that this line is representative 
of the amount of line-of-sight gaseous extinction. 

The agreement between the values obtained in the present 
work and those obtained in our previous work is very encourag¬ 
ing. Therefore, we use, hereafter, the optical depth value of the 
R(0) line observed in the spectrum of IRS 16C to correct 
the data cube for the contribution of this gaseous absorption in the 
foreground molecular clouds. 


4 THE RESIDUAL CO R(0) GASEOUS AND ^^CO 
SOLID ABSORPTIONS 

The results obtained in the previous section make the estimate of 
the continuum in the spectral range quite robust. In the following, 
we use this continuum to derive optical depth maps of the residual 
absorptions. 

The gas absorption: To correct the map shown in Fig. [5] 
for the foreground R(0) gaseous line, one has to subtract 

0.1 + 0.02. The residual optical depths imply then residual gaseous 
CO, possibly in the local medium. In Fig. |5l all pixels with green 
to red colours typically show residual ^^CO R(0) gaseous line. 
The pixels falling at the position of IRS IW and IRS 21 show 
optical depths between 0.1 and 0.17-0.18. This implies that the 
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Figure 5. Optical depth map of the gaseous R(0) line obtained from 
the corrected data cube for the foreground ice band and dust extinc¬ 
tion. Contours of the integrated M-band map are overlaid. 


corrected pixels for the foreground gaseous absorption lie in the 
interval 0 to 0.07+0.02 which is in agreement with the values 
found in Moultaka et al. (2009) of 0.04 for IRS IW and of 0.03 
for IRS 21. At the position of IRS 3, we find corrected values for 
the foreground gaseous absorption between 0 and 0.04+0.02 (in 
agreement with our previous analysis that resulted in an optical 
depth of 0.03). This implies that 10 to 30% of the absorption 

occurs in the local medium at the location of these sources. In 
general, it is not straightforward to compare our results with those 
of the high resolution data by Goto et al. (2014) who observed 
the above mentioned three sources of the Galactic Centre. Indeed, 
their spectral resolution of 3 km/s is more than 100 times larger 
than ours and the slit width of 0.2 arcsec they used on CRIRES 
instrument is three times narrower than ours. Despite that, our 
finding is in agreement with that of Goto et al. (2014) who show 
that about 20% of the CO absorption towards IRS 1 and IRS 3 
occur in gas local to the Galactic Centre. They find three main 
foreground absorption complexes due to cold and dense clouds 
in the three sipral and lateral arms at negative velocities and 0 
km/s. The residual absorptions outside these velocity ranges, in the 
positive side, were interpreted by the authors as being due to an 
extension of warm and dense clouds from the CND. In addition, 
they associate a broad trough at negative values in the spectra of 
H 3 lines with warm diffuse gas from the Central Molecular Zone 
(CMZ) as also pointed out by Oka et al. (2005) and Goto et al. 
(2008). 

The ice absorption: To derive the optical depth map of 
the residual solid-phase ^^CO absorption from the corrected data 
cube for the foreground extinction, we used the same continuum as 
previously and considered the band flux at wavelength A {hand at a) 
and the continuum flux {Icontinuum at a) at the same wavelength (i.e. 

»• I" ^is map, thc Optical depths 
are derived at pixels where the integrated M-band map has a S/N 
higher than 7. 

In Fig.[ 6 |we show the optical depth map derived from the corrected 
data cube for the line-of-sight extinction. In this figure, we find 


residuals of the 4.675yum band with optical depth values ranging 
from 0.1 to 1. The large optical depth values can be due to the 
gaseous ^^CO P(l) line that is located at the same wavelength as 
the solid ^^CO band, as one can see in the theoretical spectrum of 
Moneti et al. (2001) plotted at our spectral resolution in Fig. IA7I 
In that figure, we measure the theoretical ratio between the opti¬ 
cal depths of the ^^CO P(l) line and the ^^CO P(2) line located at 
4.68yum and find a ratio of 1.05. Thus, measuring larger values of 
this ratio in our data would imply residues of the solid band. 
This is why, we built the map of the optical depths of the P(2) 
gaseous line using the same continuum as for the previous maps. 
This map is shown in Fig.|7] Then, we derived the map of the ratio 
between the optical depths of the absorption band at A.615pm and 
that of the P(2) line. It is shown in Fig. [ 8 ] Let us call this ratio 
R. We have: 

^ _ '^Absorption af4.675juw 

”^12 CO P{2) 

(4) 

_ '^12co solid'^'^'^'^CO P(l) 

"^1200 P(2) 


To estimate the amount of residual ice, we derived the map 
of the optical depths of the solid absorption. This map is ob¬ 
tained considering the theoretical ratio of 1.05 between the optical 
depths of the P(l) and P(2) lines. Indeed, from equa¬ 
tion |4l we have: 

_ /D ^^2cOP(l)^_ 

CO solid - “ 77 ; J^^^CO P{2) 

^l2co P(2) (5) 

= (R — 1.05)ti2cc> p(2) 

The resulting map is shown in Fig.|9] 

From the values of the optical depths obtained in this map, we can 
derive the residual column densities of the CO ice absorption using 
the equation (Sandford et al. 1988): 


N{^^ CO solid) - '^^^CO solid^l^ 


( 6 ) 


In this equation, W is the FWHM of the ice absorption 

(in cm~^) and A, the absorption strength (A = 1.7vl0“^^ cm 
molecule"^). If we measure the FWHM of the band in the tem¬ 
plate spectrum we get 6 to 7 cm~^ ~ 0.016pm and if we assume the 
same value for the residuals, we can calculate the residual column 
densities Ni^^COsoUd)- From Fig.|9] we find optical depth values of 
the CO solid band ranging from 0 to ~0.6-0.7, this implies column 
densities N{^^COsoiid) of 0 to 21-24 10^^ cm"^. Typically, the 
to H 2 abundance ratios in dense molecular clouds are of the order 
of [^^C 0 ]/[H 2 ] = 8 10 “^, hence the associated N(H 2 ) range from 0 
to ~ 3 10^^ cm“^. Thus, if we consider only the molecular column 
density associated with the observed CO ice and the Ay/N(H 2 ) ra¬ 
tio of 1 to 2 10“^^ (Bohlin et al. 1978, Moneti et al. 2001), we get 
a visible residual extinction Ay of 3 to 5 mag across the region. 
This means that the residual CO ice is larger than the assumed vari¬ 
ation of the foreground extinction derived by Schodel et al. (2010) 
AAif = 0.3mag (corresponding to AAy = 0.88 mag if we assume 
the extinction law of Martin & Whittet 1990 and AAy = l.lmag 
assuming the law by Rieke & Lebofsky 1985). This result suggests 
that part of the residuals observed in our data can be due to local 
absorption. 

In the previous maps, we find important ice residuals in 
the IRS13-IRS2 region and in the IRS1W-IRS21 region. This 
implies large amounts of dust in the regions also observed in the 
mid-infrared data by Viehmann et al. (2006). This region is also 
very bright in the ^^CO R(0) optical depth map of Fig.|5] 

The high values in optical depths of the gas for the IRS13-IRS2 
sources go along with the deep residual ice absorption. This 
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Figure 6. Map of the optical depths of the 4.675//m band corrected for the 
line-of-sight solid contribution and dust extinction with contours of 
the integrated M-band map. 


Figure 8. Map of the optical depths ratio of the absorption band at A.blSjim 
over the P(2) corrected for the foreground solid band and dust 
extinction with contours of the integrated M-band map. 



Figure 7. Map of the optical depths of the P(2) line corrected for the 
foreground solid absorption and dust extinction. 


Figure 9. Map of the residual solid-phase absorption optical depths, 
corrected for the foreground solid band and dust extinction with con¬ 
tours of the integrated M-band map. 


implies that the high amount of gas and dust in this region results 
in a high shielding of the radiation field allowing for lower 
temperatures and correspondingly high ice abundances. 

These high gas-dust density regions are at the inner edge of the 
mini-cavity, hence they may be a result of the interaction of a wind 
from SgrA* with the minispiral as has also been suggested by 
Muzic et al. (2007, 2010). This is also supported by the observation 
of IR excess sources in the IRS13N complex that are indicative for 
star formation in the central parsec (Eckart et al. 2004, 2014, and 
Jalali et al. 2013). 


5 SUMMARY AND DISCUSSION 

In this paper we present a mid-infrared mapping of the central 
half parsec of the Galaxy. The data are obtained in the wavelength 
range going from 4.6 to SAjdm. Our mapping covers the region 
in which most of the continuum emission from the minispiral is 
observed in millimetre maps (see Kunneriath et al. 2012). The 
constructed field is very successful with an accuracy of the relative 
source position recovery of the order of 1 arc sec. 

We derived a first-order line-of-sight template extinction spectrum. 
The extinction is partly due to obscuration by dust and partly to 
absorption by the ice that forms on dust grains of the ISM. 
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We, then, corrected the data cube for this first approximation 
foreground extinction. 

The resulting data cube shows residuals of the 4.675yL/m absorption 
in the region suggesting that large amounts of CO ices are present 
in the central region. However, we do not exclude the possibility 
that the residuals can be due to larger variations of the foreground 
band extinction than currently admitted. 

On the other hand, from the resulting data cube, we derived optical 
depth maps of the ^^CO R(0) gaseous line and corrected it for 
the assumed foreground gaseous absorption. We found residual 
^^CO lines in the final map. Given the shape of the theoretical 
gaseous CO line complex, one may expect residuals of cold CO 
with temperatures as low as lOK. 

These results provide additional evidence to those shown by 
Moultaka et al. (2004, 2005), that very low temperatures might be 
present near the central black hole. These low temperatures may 
conflict with the physical conditions of the local environment. In¬ 
deed, the central parsec proved to be a very complex medium where 
seem to coexist the warm dust of the minispiral (~200K) and hot 
gas of the supernova remnant SgrA East, the strong winds and 
tidal forces of SgrA^, as well as high ionizing fluxes and stellar 
winds from the local early-type stars. However, within this com¬ 
plexity, one can also observe high density pockets suggested by 
the dusty filaments observed by Muzic et al. (2007, 2010) and the 
dust embedded sources that form bowshock shapes due to their in¬ 
teraction with the interstellar material (Tanner et al. 2002, 2005). 
The presence of such high-density regions preventing the ices to 
be destroyed as well as the fact that the travel time of molecular 
cloudlets in the central parsec is of the order of their lifetime of 
about 10^ -10^ yr (Van Loon & Oliveira 2003, Mellema et al. 1998) 
favours the required conditions for cold gas and ices to survive the 
passage through the central stellar cluster (i.e. as a transient phe¬ 
nomenon). 

Even though we demonstrate that low temperatures can survive in 
the local environment of the Galactic Centre, we cannot discard the 
possibility that the residual ices and cold gases that we observe in 
the present data are the signature of a more varying line-of-sight 
extinction. Such a variable line-of-sight extinction is not obvious in 
the data by Goto et al. (2014) since only two sources were observed 
through the R(0) gaseous line. 

To provide an ultimate argument for the presence of local cold gas 
and dust, a similar study has to be done at higher spectral resolu¬ 
tion. This program is in progress and will be presented in a coming 
paper. 
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APPENDIX A: SUPPLEMENTARY FIGURES 

In this Appendix, we show supplementary figures to help in under¬ 
standing the data reduction procedure. 
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Figure Al. Step 1 of data reduction procedure: example of an A and B consecutive chopped images in which one can see the positive and negative traces. 
Both images are formerly flatfielded, corrected for cosmic rays and distorsions. Vertically is the slit axis and horizontally, the dispersion axis. The traces are 
shifted relative to each other in frames A and B thanks to telescope nodding. The positive trace in A is at the same detector location as the negative one in 
frame B. Consequently, if we subtract B from A, we get the A-B image where the positive trace has twice the intensity as that in a single frame. This allows 
us to increase the S/N and to correct for sky emission lines. 



FR2_lN5exttel.fits_0 



Figure A2. Step 2 of data reduction procedure: once the A-B frames are built, spectra of the bright sources lying in the slit are extracted, calibrated for relative 
flux and corrected for telluric lines using a calibrator star of type AOV. For each extracted spectrum i, the ’shift’ and ’scale’ parameters (x/ and Si, respectively), 
obtained for an optimal correction of telluric lines, are saved for step 3. 
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standFR2_lshmlscalefinnormfull.fits_0 



Figure A3. Step 3 of data reduction procedure: in this step, the mean values of the ’shift’ and ’scale’ parameters (x/ and Si) are calculated for a single science 
frame. These values are then used to shift and scale the calibrator star spectrum (left image). The resulting spectrum is then stacked 1024 times along the slit 
axis as shown in the right image. We call this image the ’stacked standard star’ frame. 



Figure A4. Step 4 of data reduction procedure: here, the A-B frame (shown in (a) and obtained in step 1) is divided by the ’standard star frame’ (shown in (b) 
and obtained in step 3) to get a final reduced frame (shown in c). The final frame holds the complete reduced spectra of the sources lying in a single slit. For 
each slit position, all such reduced frames are added together to increase the S/N. This ended up with a single reduced frame per slit position. 
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Right Ascension 

Figure A5. Step 5 of data reduction procedure: in this final step, the reduced and summed frames obtained in step 4, corresponding to the 22 slit positions 
on sky, are positioned properly to recover the field of view. The vertical axis corresponds to the slit axis, hence to the declination, and the horizontal one, 
corresponds to the slit width position, hence to the right ascension. This results in a data cube in which each pixel of the field of view contains a reduced 
spectrum. 



Figure A6. The template extinction spectrum of the solid CO line. 
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Figure A7. Theoretical spectra of the cold and warm CO and *^CO gases as modelled by Moneti et al. (2001) at our spectral resolution. This figure is taken 
from Moultaka et al. (2009). 
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